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Recently, the Brownian dynamics of self-propelled (active) rod-like particles was explored to 
model the motion of colloidal microswimmers, catalytically-driven nanorods, and bacteria. Here, 
we generalize this description to biaxial particles with arbitrary shape and derive the corresponding 
Langevin equation for a self-propelled Brownian spinning top. The biaxial swimmer is exposed 
to a hydrodynamic Stokes friction force at low Reynolds numbers, to fluctuating random forces 
and torques as well as to an external and an internal (effective) force and torque. The latter 
quantities control its self-propulsion. Due to biaxiality and hydrodynamic translational-rotational 
coupling, the Langevin equation can only be solved numerically. In the special case of an orthotropic 
particle in the absence of external forces and torques, the noise-free (zero-temperature) trajectory 
is analytically found to be a circular helix. This trajectory is confirmed numerically to be more 
complex in the general case involving a transient irregular motion before ending up in a simple 
periodic motion. By contrast, if the external force vanishes, no transient regime is found and the 
particle moves on a superhelical trajectory. For orthotropic particles, the noise-averaged trajectory 
is a generalized concho- spiral. We furthermore study the reduction of the model to two spatial 
dimensions and classify the noise-free trajectories completely finding circles, straight lines with and 
without transients, as well as cycloids and arbitrary periodic trajectories. 
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I. INTRODUCTION 

In the traditional description of colloidal particles, 
their shape is assumed to be either spherical or rod- 
like [IH3]) or - in other words - the particles are either 
isotropic or uniaxial, i.e., rotationally symmetric around 
a figure axis [3H5J- In the isotropic case, the location of 
the particle is described by its center-of-mass position, 
while for uniaxial particles an additional unit vector is 
needed to describe its orientation. Using various prepa- 
ration techniques, by now, it is possible to prepare col- 
loidal particles with a more complex shape than spherical 
and uniaxial in a controlled way [9 - 14J. In particular, it is 
possible to synthesize particles with an orthotropic shape 
as, for example, board- like colloids [T5]. However, while 
the theory of Brownian dynamics for spherical and uni- 
axial colloidal particles is quite advanced [TJ [2] [TCI [T7] , 
much less is known for biaxial particles which need an 
additional angle to describe their location in space (on 
top of the center-of-mass position and the orientational 
unit vector). This additional degree of freedom compli- 
cates the description of Brownian motion considerably 

PUCE]. 

A second recent development in colloid science is to 
make the colloidal particles self-propelling such that 
they are moving in space on their own (so-called ac- 
tive particles) |20[ I21j . There are numerous examples 
of self-propelling colloids including catalytically driven 
nanorods [22H25] and thermo-gradient driven Janus- 
particles [26], not to speak about other realizations of 
swimmers in vibrating granulates [27], magnetic beads 
[25] , and real biological systems [2"9"r[3T] . where, for ex- 



ample, protozoa use cilia, flagella, helical traveling waves, 
and protoplasmic flow for locomotion [35]. The simplest 
description of self-propulsion is modelled by an inter- 
nal effective force which provides an constant propulsion 
mechanism on top of the Brownian motion of the parti- 
cle |39j . In two spatial dimensions, it was shown recently 
that the simultaneous action of an internal force and an 
internal torque leads to circle-swimming 27, 140j and in 
three spatial dimensions, the behavior of the mean square 
displacement was calculated analytically [TlUIS] . 

All considerations for Brownian swimmers were done 
hitherto for uniaxial particles. We are not aware of any 
study for biaxial swimmers except for a recent model- 
ing by Vogel and Stark [33]. This is important for at 
least two reasons: first, real particles are in general not 
uniaxial and therefore the effect of biaxiality needs to 
be studied. Second, it is of general importance to gen- 
eralize the equations of a "Brownian spinning top" to- 
wards self-propulsion in order to understand and predict 
its motion on a fundamental level. Surprisingly, while the 
spinning top in an external field governed by the New- 
tonian (or Eulerian) equations is a standard reference 
model in classical mechanics, there are only few studies 
for overdamped Brownian motion of a passive spinning 
top [451146]. 

In this paper, we derive the Langevin equation for 
a biaxial self-propelled particle with arbitrary shape 
that may even be screw-like, implying a translational- 
rotational coupling [37] |3H] . The biaxial swimmer is ex- 
posed to a hydrodynamic Stokes friction force at low 
Reynolds numbers, to fluctuating random forces and 
torques, as well as to an external and an internal (effec- 
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tive) force and torque. The internal forces and torques 
control the translational and angular propulsion velocity 
and are constant in the body-fixed frame while the exter- 
nal force and torque is constant in the (space-fixed) lab 
frame. 

Due to biaxiality and hydrodynamic translational- 
rotational coupling, the Langevin equation can only be 
solved numerically. In the special case of an orthotropic 
particle, which has no translational-rotational coupling, 
the noise-free (zero-temperature) trajectory is analyti- 
cally found to be a circular helix in the absence of any 
external force and torque. This trajectory is confirmed 
numerically to be much more complex in the general case. 
Typically, there is an irregular transient motion before 
the particle ends up in a simple periodic motion. Such 
helix-like trajectories are indeed typical for swimming mi- 
croorganisms [38, 49 . By contrast, if the external force 
vanishes, no transient motion shows up and the particle 
exhibits a superhelical trajectory. For orthotropic parti- 
cles, the noise-averaged trajectory is studied numerically 
and found to be a generalized concho-spiral |50j with a 
"snail-shell" structure. We furthermore study the reduc- 
tion of the model to two spatial dimensions where there is 
only one orientational angle and the internal and external 
torques can be combined to a single effective torque. In 
this simpler two-dimensional limit, we classify the noise- 
free trajectories completely. Circles are found in the ab- 
sence of the external force. In general, the trajectories 
are straight lines with and without transients as well as 
cycloids and arbitrary periodic curves. 

This paper is organized as follows: in Sec. [IT] we present 
the Langevin equation for a general self-propelled rigid 
biaxial Brownian particle in an unbounded viscous fluid 
with low Reynolds number being at rest at infinity. Sec- 
is dedicated to special analytical solutions and 
to more general numerical calculations for the 
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tion 
Sec. 

Langevin equation. Finally, we conclude and give an out- 
look in Sec. El 



II. LANGEVIN EQUATION 

In this section, we describe the Brownian motion of a 
biaxial self-propelled particle suspended in an unbounded 
viscous fluid at rest at infinity for low Reynolds numbers. 
It is assumed that the colloidal particle is rigid and has 
a constant mass density. The motion of this particle is 
characterized by the translational center-of-mass veloc- 
ity r = dr/dt with the center-of-mass position r(t) and 
the time variable t as well as by the instantaneous angu- 
lar velocity Co{t). The Brownian motion of colloidal par- 
ticles with arbitrary shape involves a coupling between 
the translational and the rotational degrees of freedom, 
which was described theoretically, for example, by Bren- 
ner [47l|48|. In 2002, Fernandes and de la Torre [IS] have 
proposed a corresponding Brownian dynamics simulation 
algorithm for the motion of a passive rigid particle of ar- 
bitrary shape. The underlying equations of motion were 



generalized to an imposed external flow field for the sur- 
rounding fluid by Makino and Doi |46) . 

Here we appropriately generalize this description to 
internal degrees of freedom, i. e., to a self-propelled biax- 
ial particle that experiences an internal effective force Fq 
and torque To, which are both constant in the body- fixed 
system. This models a biaxial microswimmer. Of course, 
a swimmer is in principle force-free and torque-free but 
the internal forces and torques are meant to be effective 
quantities which govern the propulsion mechanism of the 
particle. Using a compact notation, we introduce the ba- 
sic completely overdamped Langevin equation for three 
spatial dimensions 



D = /32?(f) (ft- 1 (f)i? - V ? U(f) + n-\i)k) 



(1) 



for a self-propelled Brownian spinning top. Let us first 
explain the notation step by step. The biaxial par- 
ticle has the position r = (xx, a?2, #3) and the orien- 
tation w = (<f),6,x)i which is given in Eulerian an- 
gles [51]. We summarize translational and rotational 
degrees of freedom by a compact 6-dimensional vector 
y = (r, vj) , which obviously involves a generalized veloc- 
ity = (r, Co) (with Co denoting the angular velocity) and 
a gradient = (Vf, V^). The latter gradient is com- 
posed both of the usual translational gradient operator 
W = {d Xll d X2 ,d X3 ) acting on the Cartesian coordinates 
of r and the rotational gradient operator = iL [52], 
given by the product of the imaginary unit i and the an- 
gular momentum operator L = (Li,L2,L3) in Eulerian 
angles. In the space-fixed coordinate system, the angu- 
lar momentum operator L is given in Eulerian angles by 



iL\ = — cos(0) cot(9)dj, — sm((/))dg 

+ cos(</>) csc(0)<9 x , 
iL 2 = — sin(0) cot(0)d^ + cos(0)<9g 

+ sin(0) csc(9)d x , 
iL 3 = d<p . (4) 



(2) 
(3) 



The angular velocity Co is expressed in Eulerian angles by 
C3 = M(m)w (5) 

with the tensor [53] 

(0 — sin(</>) cos(0) sin(0)\ 
cos(^) sin(0) sin(0) , (6) 
1 cos(0) / 

/- cos(0) cot(0) - sin(0) cot(0) l\ 
M" 1 ^) = -sin(0) cos(0) (7) 

y cos(</>) csc(0) sia((f>) csc(0) OJ 

and the Eulerian angular velocities w — dm /At. Further- 
more, in Eq. ([!]) the compact notation Kq — (Fq, Tq) for 
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the generalized force is used, which combines the effec- 
tive propulsion force and torque. In general, the biaxial 
particles is also exposed to an external potential U(f) 
giving rise to an external force -Fext = — Vj?C/ and an 
external torque T ext = — Vgf, which we both consider 
to be constant in the sequel in order to keep the model 
simple. 

The 6x6-dimensional matrix 7\L -1 (y) is associated with 
the geometric transformation from the body-fixed frame 
to the lab frame. Its inverse TZ(f) is given as the block 
diagonal matrix 



K(f) =diag(R(ro),R(w)) 
with the submatrices 

R(w) = R 3 ( X )R 2 (6)R 3 (<j>), 
R- 1 (m)=R T (m)=R 3 (- ( f ) ) R 2 (-0)R 3 (- X ) 



(8) 



(9) 



where the elementary rotation matrices Ri(<p) describe a 
clockwise rotation (when looking down the axes) around 
the ith Cartesian axis by the angle ip for i £ {1,2,3}: 



'10 
— ( cos(cp) sin(yj) 
,0 — sin(<^) cos(f) . 

^cos((p) — sin(<p) N 
R 2 {tp) =(01 

v sin(</?) cos(</?) 

cos(ip) sin(^j) s 
R 3(v) = ( -sin(p) cos(ip) 
1, 



(10) 

(11) 

(12) 



The particle shape and its hydrodynamics enter in the 
generalized short-time diffusion (or inverse friction) ten- 
sor 2?(y). This can be expressed as the 6 x 6-dimensional 
matrix 



D rr (w) D (w) 
D RT (w) D RR (w) 



1 



K- 1 ®U- l K<$) 



(13) 



where D TT (w), D TR (w) = (D RT (w)) T , and D RR (w) are 
3 x 3-dimensional submatrices that correspond to pure 
translation, translational-rotational coupling, and pure 
rotation, respectively [54]. Here, r\ is the viscosity of the 
embedding fluid and j3 — l/(k^T) with the Boltzmann 
constant k-Q and the absolute (effective) temperature T 
denotes the inverse thermal energy. The matrix % is 
constant. It only depends on the shape and the size of 
the Brownian particle and is composed of the symmetric 
translation tensor K, the not necessarily symmetric cou- 
pling tensor Cs |55) . and the symmetric rotation tensor 
n s 081156]: 



H = 



K C| 

c s n s 



The Langevin equation ([!]) thus involves altogether 21 
shape-dependent parameters. 

Finally, k(t) = (/o,to) denotes the stochastic force 
fo{t) and torque To(t) due to thermal fluctuations, that 
act on the Brownian particle in body-fixed coordinates 
for T > 0. This thermal noise k(t) is assumed to be 
Gaussian white noise with mean 



and correlation 



(*(*)) = o 



, 2 V 



(k(ti) ® k(t 2 )) = S(ti - t 2 ) 



(15) 



(16) 



(14) 



where ( • ) denotes the noise average. 

We remark that the dynamics that is described by 
the Langevin equation depends on the definition of 
the stochastic contribution. In the present form, the 
Langevin equation is only valid, if the multiplicative noise 
is defined in the Ito sense [ST]. Then, the additive drift 
term Vj • T>(f) at the end of the Langevin equation guar- 
antees that the solutions of the Langevin equation respect 
the Boltzmann distribution, when the system is in equi- 
librium at temperature T. Other definitions than that of 
Ito like, for example, the Stratonovich formulation, re- 
quire the addition of further special drift terms to the 
right-hand-side of the Langevin equation |58J. This cir- 
cumstance is always relevant in the case of multiplicative 
noise, but the necessity of the adaptation of the Langevin 
equation to the definition of the stochastic noise has been 
missed in previous work [451 146j . 



III. SPECIAL ANALYTICAL SOLUTIONS OF 
THE LANGEVIN EQUATION 

The Langevin equation ([T]) represents a system of six 
coupled nonlinear stochastic differential equations [SH] 
that cannot be solved analytically in general. There 
exist only a few analytical solutions for rather special 
situations. Several simple Langevin equations for self- 
propelled spherical or uniaxial particles in two or three 
spatial dimensions are known from literature [401442] and 
appear to be special cases of our general Langevin equa- 
tion ([lj. These Langevin equations are not as general 
and complicated as Eq. ([T]) and can be solved analyti- 
cally. Other analytically solvable special cases of Eq. 
are obtained for orthotropic particles in the absence of 
thermal fluctuations. 



A. Three spatial dimensions 

In the general three-dimensional case, it is not even 
possible to solve the Langevin equation analytically 
if the stochastic noise is neglected. Further simplifica- 
tions that reduce the number of the degrees of freedom 
or diagonalize the matrix H are necessary in order to 
obtain soluble cases. 
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B. Two spatial dimensions 



and the coupling vector 



An obvious simplification of the general Langevin 
equation ([I]) is the restriction to two spatial dimensions. 
A two-dimensional analog of the general Langevin equa- 
tion can be obtained by choosing x 3 — 0, 6 = tt/2, 
and x = 0; leaving only a single azimuthal angle <p. In 
analogy to the notation in Ref. [JDj, we further define 
K-^Ko = (F A ,0,0,0,M), k(t) = (0,/x,/||,-r,0,0), 
and Tl^iffk = (/,0,0,0,t), where F A = -F||W|] + 
F±_u±_ is the internal driving force and M is the in- 
ternal driving torque. The orientation vector un = 
(cos(</>), sin(^)) denotes the figure axis of the particle and 
u± = (— sin((/>), cos(</))) is its orthogonal complement. 
Similarly, the vector f(t) denotes the stochastic force and 
r(t) is the stochastic torque acting on the particle. The 
stochastic force vector f(t) is decomposed like the inter- 
nal driving force: f(t) = f\\u\\ Moreover, the new 
two-dimensional position vector is r = (x\, X2), the corre- 
sponding gradient is V? = (d Xl ,d X2 ), and the stochastic 
noise is characterized by = fu , — r) with mean 



and correlation 



(17) 



(18) 



where ~H = (Hij)i,j=2,3,4 is a 3x3-dimensional submatrix 
of %. The Langevin equations for two spatial dimensions 
are then given by 

f=B l + (3(D T (F A - VJJ + f) - D C (M - d^U + r)) , 
4> = P(D K (M - d^U + t)-D c - (Pa ~ V?U + /)) (19) 

with the drift vector 

=B l }u {] +Bj-u x , (20) 

which is in accordance with the interpretation of the so- 
lution of the Langevin equations as an Ito process [50] . 
the translational short-time diffusion tensor 



D T (</>) = Diu\\ <g> Uj| + D 2 (u\\ ® uj_ + uj_ ® ii||) 
+ D 3 ii ± <g> u ± , 



(21) 



(22) 



They involve only 8 instead of 21 shape-dependent pa- 
rameters. These are the translational drift coefficients 



^^(r) M -r) 15 ) 
1 



p v \v tie v- '..I 
the translational diffusion coefficients 



1 



(«-') 



the coupling coefficients 



1 



33 



1 

fa 



^12 ' 



^(^ ^23 
1 



(u- 



1 13 



and the rotational diffusion coefficient 



'33 



(23) 
(24) 

(25) 
(26) 
(27) 

(28) 
(29) 

(30) 



Some of these eight coefficients are zero or equal, respec- 
tively, if the described Brownian particle is symmetric. 
Table [T] gives an overview about possible symmetries of 
the particle's shape and the corresponding properties of 
the shape-dependent coefficients ( p3] ) -(|30| ). 

Although the Langevin equations ( |19[ ) for two spatial 
dimensions are simpler than Eq. 0, they are still cou- 
pled nonlinear stochastic differential equations and thus 
not analytically solvable. However, if the external po- 
tential U(r, 4>) is set to zero, the Langevin equations 
can be solved analytically and the center-of-mass tra- 
jectory becomes a circle (T = 0) or a logarithmic spiral 
(T > 0) like in Ref. g0]. The analytical solution for 
U(r, 0) = f\\ — fj_ — t — is given by 



f(t) =f + ^(F\\D 2 + F ± D 3 - MD^)(u\\ (fa + ut) - u\\ (fa)) 
— (F\\Di + F ± D 2 - MDl)(u ± (<p +ut) - u±(fa)) 



LO 

<f>(t) = fa+Ljt 



(31) 
(32) 



with the angular velocity 

w = /3(MD R - D^F {1 - D^F ± ) 



the initial position fo = r(0), and the initial orienta- 



(33) 
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TABLE I. Connection between the symmetry of the particle shape and the parameters ( 23 (- ( 30 ( in the Langevin equations 
( 19 1 for two spatial dimensions. 



type of shape: 


uniaxial 


uniaxial uniaxial 


uniaxial 


isotropic 


symmetries: 


no symmetry 


inflection symmetry inflection symmetry 


inflection symmetry 


rotational symmetry 


invariance 
properties: 




Xl — > —Xl X2 — > —X2 


xi — >• —si, 
£2 — > —2:2 


-> + A0 
V A(f> 6 [0, 2tt) 


shape-dependent 
parameters: 

Di / 0, Dr. / 


b] 1 ± 0, Bt ± 0, 
D 2 ^ 0, D 3 ^ 0, 
Dl ? 0, * 


b\ = 0, Bi ± 0, b|' £ 0, Bf 1 = 0, 
D 2 = 0, D 3 / 0, D 2 = 0,D 3 =£ 0, 
Z>Uj / 0, Z?£ = Z)|L = 0, Z>£ ?4 


bJ 1 = 0, Bt = 0, 
B 2 = 0, As ^ 0, 
DjL = 0, £>£ = 


bJ 1 = 0, Bf 1 = 0, 
D 2 =0, B 3 = Di, 
b| = 0, Dh = 



tion <fio = (f>(0). If instead of U(f,(j)) the coupling coef- 
ficients Dq and vanish, the Langevin equations be- 
come similar to the Langevin equations for the Brown- 
ian circle swimmer in Ref. |40j . but with an additional 
non-diagonal coefficient D2 in the parametrization of the 
translational diffusion tensor Dt(</>) and a more general 
driving force which is not necessarily parallel to the 
figure axis. For D2 = = Dq. = F± = 0, the Langevin 
equations ( 19 1 are equivalent to the Langevin equation in 
Ref. |40j . An additional constraint on spherical particles 
that are only able to move along the x\ axis leads to the 
Langevin equations for a spherical self-propelled particle 
on a substrate [41] . With similar simplifications it is also 
possible to derive the Langevin equations for spherical or 
anisotropic uniaxial self-propelled particles in two spatial 
dimensions that are discussed in Ref. 42 . 



tions for orthotropic particles: 

?= /3D TT (ro)(R- 1 (ro) J F - V P U + R-\m)k 1 ) , 
Q = /3D RR (w)(R- 1 (w)f - V*t7 + R _1 (*)jfe a ) (34) 
+ V,s • D RR (n7) . 

Here, the Gaussian white noises k\ and k 2 are inde- 
pendent and defined as the first and second part of 
k = (fcijfca), respectively. The Langevin equations for 
spherical or uniaxial particles, that are considered in Ref. 
[42] , are special cases of the more general Langevin equa- 
tions (34) for biaxial orthotropic particles. To be able to 



solve these Langevin equations analytically, it is at first 
necessary to neglect ki and k 2 , i.e., to consider the case 
T = 0. A further negligence of the drive or the external 
potential leads to two special cases which are analytically 
solvable. 



C. Orthotropic particles 



Settling orthotropic passive particle: 



Another possibility to simplify the Langevin equa- 
tion ([!]) considerably is the exclusive consideration of 
orthotropic particles. All geometric bodies with three 
pairwise orthogonal planes of symmetry like spheres, 
spheroids (ellipsoids of revolution), biaxial (or triaxial) 
ellipsoids, cylinders, cuboids, and some prisms belong 
to this important class. In general, orthotropic parti- 
cles have no translational-rotational coupling so that the 
coupling tensor Cs vanishes. Furthermore, the trans- 
lation tensor K and the rotation tensor Sis are diag- 
onal for orthotropic particles. These properties of K, 
Cs, and Sis c& n be derived from the circumstance that 
the center of mass, which should be chosen as reference 
point S, and the mutual point of intersection of the three 
planes of symmetry of an orthotropic particle coincide 
[55] , The vanishing of the coupling tensor Cs results 
from the fact, that the point of intersection of the three 
pairwise perpendicular planes of symmetry of the parti- 
cle is identical with the center of hydrodynamic reaction 
for orthotropic bodies. With these considerations, the 
Langevin equation ([I]) simplifies to the Langevin equa- 



A particle without drive, i.e., with Fq — Tq = 0, 
moves only under the influence of the external poten- 
tial U(f, zu). In the case of a constant gravitational field, 
only the constant external force F oxt = —VpU acts on the 
particle and the external torque — V^t/ vanishes. The 
motion of such a sedimcnting particle is well known from 
literature [56]. It is characterized by a constant velocity 
r and a constant orientation id: 



f — /3D TT (roo)-F e xt = const, 
w = TJ7o = const. 



(35) 



2. Self-propelled orthotropic particle: 

If the external potential U (r, w) is neglected instead 
of the drive in Eqs. (34) for T = 0, they describe the 



helical motion of an arbitrary orthotropic self-propelled 
particle in the absence of external and random forces and 



torques: 



V 

LU = R -1 ^)^ 1 -^ . 

7] 



(36) 



These equations of motion are trivial in the body-fixed 
coordinate system, where the velocities r and lu are con- 
stant. Since the angular velocity lu with the initial orien- 
tation ro = ra(0) is constant in body-fixed coordinates, 
it is also constant in the space-fixed system. This means 
that one expects a helix for the center-of-mass trajectory 
as it is known from the motion of protozoa like Euglena 
gracilis [38] and bacteria like Thiovulum majus |49j . The 
analytical solution of Eq. (p6|) is in fact the circular helix 



(37) 



-./ - (Co X Vq) X LU . .. _, . 
r(t) = ro + ~ sin(|w|t) 



LU X Vq 



UJ 



with axis 



along LU, radius 



and pitch 



= r + 



i(\Lu\t)) + —^LUt 



LU X Vq 



LU 



KK- 1 ^) x (figif )| 



?7 2 1 c<3 1 2 



h = 2tt 



\(K- 1 F ) ■ (fig 1 T )| 



?7 2 |u;| 2 



(38) 



(39) 



(40) 



where r*o = r(0) is the initial position, vq — r(0) = 
R~ 1 (zuq) K _1 Fo/r/ is the initial velocity, and |u;| = 
IJlg 1 ^/^! is the modulus of the angular velocity. This 
helical trajectory is shown schematically in Fig.[T] When 
also a constant gravitational field is taken into account, 
the helix becomes deformed and, for example, its cross 
section might become elliptic, but the axis of the helix 
remains a straight line. In the case T > 0, where we have 
to consider the stochastic contributions in the Langevin 
equations (34), it is, however, no longer possible to find 



analytical solutions. This case requires the usage of ap- 
propriate numerical integrators for stochastic differential 
equations and is studied in the next paragraph. 



IV. NUMERICAL CALCULATIONS 

In general situations, where analytical solutions do not 
exist, the Langevin equation ([!]) can only be investi- 
gated with the help of numerical methods. Appropri- 
ate numerical methods in increasing order of the trun- 
cation error are the Euler-Maruyama method, the Mil- 
stein method, and stochastic Runge-Kutta methods [50] . 




FIG. 1. The center-of-mass trajectory for T = and a con- 
stant external potential is a circular helix with radius r and 
pitch h that evolves from the rotation of the orthotropic par- 
ticle with the constant angular velocity w. 



which have to be derived in the Ito sense in order to be 
compatible with the Langevin equations in this paper. If 
there are no thermal fluctuations (T = 0), the stochas- 
tic Langevin equation ([!]) becomes deterministic and a 
standard Runge-Kutta scheme of high order can be ap- 
plied. The numerical results for T = and T > 0, that 
are presented in the following, have been obtained by an 
explicit fourth-order deterministic Runge-Kutta scheme 
|61H64j and by a multi-dimensional explicit stochastic 
Runge-Kutta scheme of weak order 2.0 for Ito stochas- 
tic differential equations |60j . respectively. During the 
whole section, Fq, To, F cxt = — Vpf/, and T ext = — V^U 
are assumed to be constant vectors, which do not de- 
pend on r or w. Furthermore, arbitrary Brownian par- 
ticles with a hydrodynamic translational-rotational cou- 
pling and orthotropic particles without a translational- 
rotational coupling in two and three spatial dimensions 
are considered for T = and T > 0. In parallel to the 
previous section, this section is divided into a first sub- 
section about the general Langevin equation for three 
spatial dimensions, a second subsection about two spa- 
tial dimensions, and a third subsection about orthotropic 
particles. 



A. Three spatial dimensions 

For arbitrarily shaped particles in three spatial dimen- 
sions, one finds various differently shaped trajectories as 
solutions of the Langevin equation (JlJ . Figure 2] gives 
a selection of typical trajectories that can be observed 
for arbitrarily shaped particles with an arbitrary drive at 
T = 0. In order to sample some typical solutions, we have 
randomly chosen the 21 shape-dependent parameters for 
the matrix "H, the components of the internal force Fq 
and torque To, the components of the external force F ext 
and torque T ext , and the initial conditions. More than 
100 random parameter combinations were considered. In 
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. transient regime . 




(a) Fext / 0, Text / 



(b) Fext / 0, Text = 





no transient regime 

(c) Fext = 0, Text / (d) Fext = 0, f ex t = 

FIG. 2. Typical trajectories of arbitrarily shaped active par- 
ticles in three spatial dimensions for constant vectors F$ 7^ 0, 
To 7^ and temperature T — 0. The external force T ex t = 
— Vft/ and torque T cxt = — V#t/ are constant, too. For a 
non-vanishing external force, the particle's center-of-mass tra- 
jectory starts with an irregular transient regime and changes 
into a periodic motion, as it is shown in plots (a) and (b). 
The general periodic motion (a), that is observed if there is 
also a non-vanishing external torque, reduces to a circular 
helix (b) parallel to the external force, if there is no external 
torque. The other two plots (c) and (d) show the situation for 
a vanishing external force, where a transient regime is not ob- 
served. There, the trajectory is either a superhelix-like curve 
(c) parallel to a non-vanishing external torque or a circular 
helix (d), if there is no external torque. However, the tra- 
jectories (a)-(c) can also be irregular. Straight trajectories, 
that are preceded by a transient regime if T cxt 7^ 0, can be 
observed, too. 



doing so, four different cases with vanishing and non- 
vanishingvectors for _F ext and T ext were distinguished 
(see Fig. |2j). Depending on the choice of -F cxt an d ^ext, 
the observed trajectories appeared to share common fea- 
tures and to be distinguishable into four different classes. 
The four trajectories that are shown in Fig. [2] are repre- 
sentatives of these classes. They can be characterized as 
follows: If there are both an external force and an exter- 
nal torque, self-propelled particles that start with an ir- 
regular transient regime and end up in a simple periodic 
center-of-mass trajectory are usually observed [see Fig. 
[2]ja)]. Note that the length of the initial transient regime 
as well as the periodicity length of the final periodic mo- 



tion depend on the particular parameters and may be- 
come rather long. In the case of no external torque, the 
periodic motion after the transient regime is a circular 
helix with its axis being parallel to the direction of the 
external force vector. This situation is illustrated in Fig. 
|2j;b). The analog case of an external torque but no ex- 
ternal force is schematically shown in Fig. [2jc) . There, a 
superhelix-like curve with the orientation parallel to the 
direction of the external torque vector and without a pre- 
ceding transient regime is observed. In contrast to the 
trajectory in Fig. [2]ja), the complicated superhelix-like 
curve does not turn into a simpler periodic curve after 
some time, since there is no transient regime for F cxt = 0. 
As fourth case, the motion in the absence of both exter- 
nal forces and torques is shown in Fig. [2^d) . It appears to 
be a circular helix. Also in this case, a transient regime is 
not observed. In the situation of Fig.[2]Ja) and (b), com- 
pletely irregular trajectories appear, when the transient 
regime is very long. Even in the transient-free situation 
of Fig.[2][c), irregular trajectories are observed, when the 
rotational frequency ratio between the immanent rota- 
tion of the self-propelled particle and the rotation due to 
the external torque is irrational. Furthermore, straight 
trajectories appear as a special case and are preceded by 
a transient regime if F cxt 7^ 0. A complete and detailed 
classification of all trajectories that can be observed in 
three spatial dimensions is, however, not possible, since 
the number of the parameters that define the shape of 
the particle and all further relevant quantities like inter- 
nal and external forces and torques is quite big, but this 
number is much smaller in two spatial dimensions, where 
a more detailed classification is possible. 



B. Two spatial dimensions 



The Langevin equations ( 19 ) for two spatial dimcn 



sions were solved analytically in Sec. IIIB for U(f, <f>) = 
and T = 0. Here, we consider the case T = 0, too, 
but now for a constant non-vanishing external force 
-Fcxt = — V^J7, since the external force leads to various 
different non-trivial trajectories. The observed trajecto- 
ries are classified with respect to the shape and the kind 
of self-propulsion of the particle in Tab.|n] Since the drift 
coefficients Bj and B^~ can be neglected for T = 0, the 
shape of the particle is only described with the six param- 
eters D±, Z?2, D3, Dq, Dq, and T>r, where D\ and £>r 
can be set to one by a suitable rescaling of the length and 
time scales. For the remaining parameters, particles with 
a translational-rotational coupling and particles without 
a translational-rotational coupling as well as asymmetric 
particles, particles with one axis of symmetry, particles 
with two mutually perpendicular axes of symmetry, and 
isotropic particles with rotational symmetry are distin- 
guished. Moreover, the constant parameters F±, and 
M e g — M — dff,U are used to describe the self-propulsion 
of the particle. Four situations of a non-vanishing inter- 
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nal force Fa and torque M, a drive only by either an 
internal force or an internal torque, and a passive par- 
ticle with a vanishing drive are considered. A further 
distinction with respect to the external force and torque 
is not necessary, because the external force can always be 
chosen bigger than zero, since the case of a vanishing ex- 
ternal force has been proven to lead to a trivial circular 

-d A U 



trajectory in Sec. 



IIIB 



and the external torque 
is already included in the effective torque M e g. In gen- 
eral, straight lines with an aperiodic transient regime, 
arbitrary periodic curves, cycloids, and simple straight 
lines were found as trajectories in two spatial dimensions 
by random choices of the parameters. These trajectories 
have still the basic features of their three-dimensional 
analogs, but are much simpler to describe. It is only for 
particles with translational- rotational coupling, i. e., par- 
ticles where at least one of the coefficients Dq and Dq 
does not vanish, that the straight trajectories are pre- 
ceded by an initial transient regime. These trajectories 
are characterized by a monotonous rotation of the par- 
ticle when it starts moving and an ensuing rotation-free 
straight motion. In the transient regime, the particle ro- 
tates until the internal torque, the external torque, and 
the additional torque due to the translational-rotational 
coupling compensate each other. For active particles 
with translational-rotational coupling, also periodic tra- 
jectories are observed, where a canceling of the total 
torque does not happen during the initial rotation. This 
is not the case for symmetric particles with a vanishing 
effective torque M e g and for passive particles, for which a 
periodic trajectory is not observed. The motion of sym- 
metric particles without translational-rotational coupling 
is always periodic or constant. In both cases, the trajec- 
tory is parallel to the direction of the external force, if 
the effective torque M c g is not zero. Without the ef- 
fective torque, only straight trajectories are observed for 
these particles. Solely in the case of passive rotation- 
ally symmetric particles, the orientation of these straight 
trajectories is parallel to the external force. 



C. Orthotropic particles 

To regard the influence of thermal fluctuations on 
the motion of a biaxial self-propelled Brownian parti- 
cle, the Langevin equations (34) for orthotropic particles 

For 



are solved numerically for T > in this section, 
this purpose, first of all, characteristic quantities for the 
length scale, time scale, and force scale are chosen and 
the Langevin equations ( 34 ) are rescaled to dimension- 



less units. A suitable choice for the characteristic length 
l c , the characteristic time t c , and the characteristic force 
F r is 



/ A max (D TT ) 

A m ax(D RR ) ' 
1 



A max (D RR ) 



(41) 
(42) 



tc. 



»?A max (D li ) , 



(43) 



where A max ( • ) denotes the biggest eigenvalue of the 
respective matrix. These characteristic quantities are 
used to express the position f = r*7 c , time t = t't c , 
forces P j= F^F C , F cxt = F e ' xt F c , torques f = f{,FJ c , 
T cx t — F ' xt F c / c , translation tensor K = K'l c , and ro- 
tation tensor f2s = ^s^c by the dimensionless position 
f = (x' 1 ,x' 2 ,x' 3 ), time t', forces Fq, F c ' xt , torques Tq, 
Tg Xt , translation tensor K', and rotation tensor fig, re- 
spectively. In the rescaled Langevin equations, the pa- 
rameter 



T' 



2t c 




; ( D RR) 



C (D TT ) 



oc T 



(44) 



appears as a dimensionless temperature. This param- 
eter is varied and fluctuation-averaged trajectories are 
calculated for different temperatures with fixed initial 
conditions fg and w' . The results for the case of van- 
ishing external forces and torques, where the trajectory 
for T' = T = is known from the analytical solution 
in Sec. |III C| to be a circular helix, are shown in Fig. 
[3] For this figure, we chose the dimensionless forces 
F^ = (-0.5,0,3), F c ' xt = 0, torques = (-1,0,0), 



0, translation tensor K' = diag(l,2,3), rota- 



tion tensor Q' s = diag(l,3,4), initial conditions Tq — 

04,0) ^2,0 j^io) = 0, w' = w = (0,7r/2,0), and tem- 
peratures T' — 0,0.05,0.1,0.3. It is apparent that the 
center-of-mass trajectory for T' = and the fluctuation- 
averaged center-of-mass trajectories for T' > have no 
transient regime in Fig. [3] This is also the case for 
non-vanishing constant external forces and torques and 
a general feature of orthotropic particles in contrast to 
less symmetric particles with a translational-rotational 
coupling. In the presence of thermal fluctuations, the 
helical motion of the self-propelled orthotropic particle 
is damped exponentially with time and the fluctuation- 
averaged center-of-mass trajectory becomes a concho- 
spiral |50j . whose radius and pitch decay exponentially 
with time [see Fig. |4^a) and (b)] . This result was con- 
firmed by a fit of the numerical solutions with the gen- 
eral parametrization of a concho-spiral and agrees with 
the observation of a logarithmic spiral, also named spira 
mirabilis by Jacob Bernoulli, in the two-dimensional spe- 
cial case of our Langevin equations, that is investigated 
in Ref. [3D]. In the situation of Fig. [3j the axes of 
the concho-spirals are parallel to the X3-axis and can be 
parametrized by [3DJ ISP] 



x[(t' 
x' 2 (t> 

4(f 



= x' 10 + a[ ( cos(0 o ) - cos(0(t'))e 7;t ' 
- ol 2 (sin(0 o )- sinOK^e-^') , 

= + "i ( sin(^o) - sH^t')^-^ 1 ') 
+ al,(cos((/> ) -cos((/)(i'))e" 7; *') , 

= 4,0+^max(l-(l-£ / *')e- 7 =*') , 



(45) 

(46) 
(47) 
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(a) T' = (b) T' = 0.05 (c) T' = 0.1 (d) T' = 0.3 

FIG. 3. Mean trajectories of a self-propelled orthotropic particle in the absence of external forces and torques for the dimen- 
sionless temperatures T' = 0, T' = 0.05, T' = 0.1, and T' = 0.3 (from left to right). In the left two plots, the trajectories are 
shown for the time interval ^ t' ^ 40, while the right two plots show mean trajectories with ^ t' ^ 80. 



<f>(t') = (f> Q + u't' 



(48) 



SCO 



with the angular frequency lj = \lj 
Eqs. (36 1] and the dimensionless fit parameters a' 11 a' 2 , 
and e 



7l) 72: 



max ) 



Equations (45), (46), and (48) 



describe a logarithmic spiral, which is the trajectory of 
the two-dimensional circle swimmer in Rcf. 1401 , while the 



parametrization (47) of the third spatial variable x' 3 (if) is 
here more general than in Ref . [50] . In Eq. ( 47 ) there is an 



additional term oc e', which makes sure that a helix is ob- 
tained as special case of the concho-spiral for T' — 0, i. e., 
for 7i = 72 = 0. This is, however, not the case for the 
parametrization in Ref. |50j . Based on the parametriza- 
tion (45 )-(|48|), radius and pitch of the concho-spirals can 



be derived. The fit parameters a[, a' 2 , and "/[ determine 
the dimensionless radius r'(t') = r(t)/l c of the concho- 
spirals: 



r'(t') =r , e-^ t ' 



r'(0) 



,/2 



(49) 
(50) 



Their dimensionless pitch h'(t') = h(t)/l c depends on the 



remaining fit parameters 7 2 , H' max , and e': 



h'{t') = H' m 



=-^*'((l-£Y)(l 

,2tt _2„ 

e — e ^ 



H 2) 



(51) 



Numerical values for radius and pitch are shown in Fig. 
[4^a) and (b). Furthermore, the height H max (T) of the 
concho-spiral and its dimensionless analog H' max (T') = 
H max (T)/l c , defined as the distance from the initial po- 
sition of the particle to its final position for t' oo 
measured along the axis of the concho-spiral, are finite 



and decrease monotonously, when the temperature is in- 
creased. For the numerical calculations that correspond 
to the results shown in Fig. [3j the inverse power law 



tfmaxCn « 0.82 T'" 



(52) 



was determined [see Fig. [4](c)] . When there is addition- 
ally a constant external force, the helix for T' = as well 
as the concho-spirals for T' > are deformed and their 
cross-sections can become elliptic. 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have studied the Langevin equation 
governing the motion of a driven Brownian spinning top 
describing a self-propelled biaxial colloidal particle. The 
particle is driven both by internal and external forces 
and torques, which are constant in the body frame and 
in the lab frame, respectively. This equation is non- 
trivial due to the geometric biaxiality and the hydro- 
dynamic translational-rotational coupling of the parti- 
cle and can therefore only be solved numerically. In the 
special case of an orthotropic particle in the absence of 
external forces and torques, the noise-free trajectory is 
analytically found to be a circular helix and the noise- 
averaged trajectory is a generalized concho-spiral. The 
noise-free trajectory is confirmed numerically to be more 
complex for translational-rotational coupling involving a 
transient irregular motion before ending up in simple pe- 
riodic motion. By contrast, if the external force van- 
ishes, no transient is found and the particle moves on a 
superhelical trajectory. We furthermore studied in de- 
tail the much simpler reduction of the model to two spa- 
tial dimensions. In two spatial dimensions, the noise-free 
trajectories are classified completely and circles, straight 
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(a) radius r'(t') = 


r(t)/l c 




. . . numerical data 




— non-standard fit ' 


H' m ^ = 16.57 




y\ = 0.06 




e' = -0.03 
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(b) pitch h'(t') = 


h(t)/l c 



+ + + numerical data 
0.82/r' fit 




(c) height H„ 



FIG. 4. (a) Radius and (b) pitch of a trajectory, that is shaped 
as a concho-spiral, decay exponentially with time t' . In the 
linear- logarithmic plots (a) and (b), the exponential decay is 
obvious for the radius, but not for the pitch, where the ex- 
ponential function has a linear time-dependent prefactor [see 
Eq. (51 1]. The numerical data (red dots following a curved 



line) for these plots were taken from the concho-spiral in Fig. 
[3^b) . (c) The height of the concho-spirals is inversely propor- 
tional to the temperature T' . The parameters for plot (c) are 
the same as in Fig. [3j but with more and different values for 
T' . In each plot, a straight blue line corresponding to Eqs. 
d49l-(|52l was fitted to the numerical data. 



lines with and without transients, as well as cycloids and 
arbitrary periodic trajectories are found. 

The Langevin equation derived here can be used as a 
starting point to describe more complex situations. As 
examples, we mention a particle in confinement of lin- 
ear channels [H I3HSS] or cylindrical tubes [271155], A 
sliding motion can be expected similar as in the two- 
dimensional reduction of our model |40) . In confinement, 
details of the propulsion mechanism are getting relevant 
since they result in different hydrodynamic interactions 
of the particle with the system boundaries [57] , 



A next level of complexity is given by an ensemble of 
swimmers, i. e., a finite number density, which can inter- 
act either directly by excluded volume or via hydrody- 
namic interactions. Recently, a dynamical density func- 
tional theory has been proposed for biaxial particles |68j 
following the lines given for uniaxial particles 69 . The 
collective properties of self-propelled biaxial particles are 
assumed to be rich including turbulent states, swarming, 
and jamming [70j 171] 



Finally, one can impose a non-vanishing solvent flow, 
like Couette shear flow |72j . and study a self-propelled 
particle there 73J. Our Langevin equation can be 
straightforwardly generalized to this situation [46] and 
may lead to new shear-induced tumbling phenomena of 
self-propelled biaxial particles. 
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TABLE II. Detailed classification of the trajectories of arbitrarily shaped particles in two spatial dimensions for T — with 
respect to the symmetries that are summarized in Tab.jl] For the external force, F* oxt 7^ is chosen, since all trajectories become 
circles otherwise. In the plots below, _F cxt is always oriented downwards, i. e., in the negative ^-direction. Internal and external 
torques are combined to the effective torque M e g = M — d$U. 



particle shape 


no symmetry 


a;i-axis or a;2-axis is 
axis of symmetry 


Xi-axis and X2-axis are 
axes of symmetry 


rotational symmetry 


D 2 ¥= 0, D 3 ± 0, 
d!1 ^ 


D 2 = 0, D 3 / 0, 
d|L / v Djk 


D 2 = 0, D 3 £ 0, 
dJJ, = 0, = 


D 2 = 0, D 3 = Di, 
Dl = 0, = 




Ik 



Tk 

6? 


Ik 




/ \ 

u 








straight line after transient 
regime or periodic curve 6 


straight line after transient 
regime" or periodic curve 6 


periodic curve c || F ex t 


cycloid c || F ext 




11 
« 

HI 

^ H 


ik 




\ 




\ 


straight line after transient 
regime" or periodic curve 6 


straight line after 
transient regime" 


straight line"" 


straight line d 




Ik 

1 



II 



11 












straight line after transient 
regime" or periodic curve 6 


str. line after trans, regime" 
or periodic curve 6 || F ex t 


cvcloid c II F ov + 

II e x t 


straight line c || F ex t 




11 
its 



II 



II 
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J 








straight line after 
transient regime 


straight line after 
transient regime" || F ex ± 


straight line" 1 


straight line d || F ex t 



a The particle rotates monotonously until it reaches its final orientation. Then the angle (j> remains constant. 
6 The angle (f> and the center of mass of the particle describe periodic curves with the same periodicity. 
c The particle rotates with a constant angular velocity, i. e., (j> tx t. 
d The orientation of the particle is constant: <f> — const. 



